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Introduction 


There are many systems which are described by nonlinear operators. A quadratic stochastic operator 
(QSO) is one of the simplest nonlinear cases. A QSO has meaning of a population evolution operator 
and it was first introduced by Bernstein in [1]. For more than 80 years, the theory of QSOs has been 
developed and many papers were published (see e.g. [4]-[7]). In recent years it has again become of 
interest in connection with its numerous applications in many branches of mathematics, biology and 
physics. 


Let = {1, 2 is m} be a finite set and the set of all probability distribution on E 


SS ix = (Xx,,%),...,x, )€0":x, 20, for any i and ) x, = i (1) 

i=l 
be the (m — 1) -dimensional simplex. A QSO is a mapping defined as V:S”"'—>S”"" of the 
simplex into itself, of the form V (x) =x’ e¢ §”", where 


5 DP x, k € E, (2) 


i,j=l 


and the coefficients P. , satisfy 


=P... 20, » P,,, =1 forall i, j ¢ E. (3) 
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The trajectory (orbit) a , of V for an initial value x) e S”” is defined by 


n>=0 
x" =y (x'”) -yo) ie), n=0,1,2,... . 


One of the main problems in mathematical biology is to study the asymptotic behavior of the 
trajectories. This problem was solved completely for the Volterra QSO. 
The operator V is called Volterra QSO, if P,, =0 for any k ¢ {i, i}, i, j,k € E. For the Volterra 


QSO the general formula was given in [4], 


KX [ + Ya, } (4) 
i=l 


where a, =2P,, —1 fori#k and a, =(0. Moreover, a,, =—a,, and \a,,| <1l forall i,k ek. 


In [4], the theory of Volterra QSO was developed using theory of the Lyapunov functions and 
tournaments. But non-Volterra QSOs were not completely studied. Because, there is no general theory 
that can be applied for study of non-Volterra operators. 


In [7] Separable Quadratic Stochastic Operators (SQSOs) where introduced. Volterra QSO (4) has a 
form as SQSO, but in [8] it was proved that it coincides with a SQSO if and only if it is a linear 
operator. 


In recent years, Cubic Stochastic Operators (CSOs) have begun to be studied, which different from 
quadratic operators. 


In this paper we consider another class of cubic operators which we call separable cubic stochastic 
operators. 


In Section 2, we recall the definition of CSOs and definitions and known results. In Section 3 for a 
SCSO defined on the two-dimensional simplex, we prove that it has three fixed points and we find 
conditions on parameter under which a fixed point is a repelling, attracting, or saddle point, the 


2” » . é 
boundary OS™ is an invariant set. 
Preliminaries and known results 


Separable quadratic stochastic operator. Let us recall some necessary definition and notations. In [7] 
Separable Quadratic Stochastic Operators were introduced as follows: The QSO (2), (3) with 
additional condition 


P= a,,D i, for all i, j,k € E (5) 


ijk 


where Wis0, €[] entries of matrices A = (a, ) and B= (,.) such that the conditions (3) are 


satisfied for the coefficients (5). 


Then the QSO V corresponding to the coefficients (5) has the form. 


x, _ (V(x)), = (A(x)), (B(x)), » (6) 
where (A(x)), = De, ; (B(x)), = Sb, X, : 
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Definition 1. [7] The QSO (6) is called separable quadratic stochastic operator (SQSO). 


From the conditions ean >0 and >, P. .= ] for all 7 : i it follows that the condition on matrices A 
k=1 

and B that a,D,, = 0, AB’ = 1, where B’ is the transpose of B and 1 is the matrix with all 

(i) 


entries 1’s. If @ 


um 


) is the 1 -th row of the matrix A and pb) = (Bread 


(0 ni5G, * is the 


J -th row of the matrix B, then from AB "=1we get 


ap? = 1), for all is J = 1,...,™m. 


For a fixed 7 , the above condition implies that 
aT 
A(b) = (1,L...1).@ 


If det (A) # (), then (7) gives pb) — p” for all k,j € E, ie., all rows of B are the same, 
therefore det (B) = (). Similarly, if det (B ) # O, then all the rows of A must be the same, so 
det (A) =0. 


Cubic stochastic operator. The CSO is a mapping W : S mY > §”" of the form 


m 


x= >? bie ae Le E, (8) 


ijk I” i- jk? 
i,j,k=l 


where Ba , are coefficients of heredity such that 


| AM ci gas pW) re i! ew 


ijkl kij.l ikj Jl kji,l jik.l j ijk Jl 


=1Vi, j,k € E.(9) 


and we suppose that the coefficients a do not change for any permutation of 1, J, k. 


Note that W is a non-linear operator, and its dimension increases with 7. Higher-dimensional 
dynamical systems are essential, but only relatively few dynamical systems have yet been analysed. 


For a given x”) eS ae the trajectory {x} of initial point x) under action of CSO (8) is 
n20 


defined by xl") = W(x" ) , where 2 =0,1,2,... with X = xO), Denote by o(x"”) the set of 


limit points of the trajectory {x} . Since {x} cS” and S"" isa compact set, it 1s 
n=0 n=0 


follows that ox) #O. If o(x"”) consists of a single point, they the trajectory converges and 
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o(x'”) is a fixed point of the operator W. A point x E S$ ” is called a fixed of the W if 


WwW (x) =X. Denote by Fix (W) the set of all fixed points of the operator W , i.e. 
Fix(W) = {x eS"! :W(x)= x}. 
Let DW (x") = (ow, / Ox, )(x’) be a Jacobian of W at the point x. 


Definition 2 ([3]): A fixed point X is called hyperbolic if its Jacobian DW (x’) has no 


eigenvalues on the unit circle in 
Definition 3 ((3]): A hyperbolic fixed point X° is called: 


(i) attracting if all the eigenvalues of the Jacobian DW (x" ) are in the unit disk; 


(ii) repelling if all the eigenvalues of the Jacobian DW (x" ) are outside the closed unit disk; 


(iii) a saddle otherwise; 
Main result 


In this section we consider CSO (8), (9) with additional condition 


P., =4a,b,c,,, forall i, j,k,LeE FE, 1) 


ijkl jl- kl? 


where Wh sDasGe E entries of matrices A = (a, ). B= (b, ) and C= (c,, ) such that the 


conditions (9) are satisfied for the coefficients (11). 


Then the CSO W corresponding to the coefficients (11) has the form 


=(W(x)), =(4()), (B09), (CC), 2 
where (A(x)), = Ya, (B(x)) = Deer and (C(x)) = Yc, 


Definition 5. The CSO (12) is called separable cubic stochastic operator (SCSO). 


Lemma 1. The condition (11) is sufficient for a CSO to be product of three linear operators, but the 
condition is not necessary. 


The CSO with coefficients P,, can be written as the product of three matrices A = (a, ). 


ikl 


B= (b,.) and C = (c,,) if and only if, for any Ly Joke € E it holds 


Les + lem at Ls ot ie + gee a ane = a,b Cy - D.C 4 aa GDC a i ‘ 
us the 


G0 jC, + 00,6, 7G, DC (13) 


jlvil pil Kl il 
condition (11) is a particular case of (13). 
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Let us consider the following matrices: 


100 - + = } fen 4 
A=|0 1 O|},B=|1l-a 1 L ice! J 1 l+a |,(18) 
0 0 1 1 l-a 1 12 1 1 

2 


where a €[—1,]1]. 


Then corresponding SCSO W : S* > S’ is: 


c= (2% +(1-a)zx, alls +X, +(1-$}s5} 


W: x =x, (x, +x, +(1-a)x,)((I+a)x, +x, +x), (19) 


ee ic ae +X, + « |(s +(1+a)x, + * 


Using the equation X, + X, + X, = 1 we rewrite the operator (19) as follows 


x = (I-ax,)(1-$s5} 


W :5 x, =x,(1-ax, )(1+ax,), @9 


a 
x, [14s Jira) 
Evidently, that if @ =O the SCSO (20) is the identity map. For this in the below, we consider the 


case when a + OQ). 


Let a face of the simplex S” be the set Tr, = {x E S -x, = 0, L~ZQAaAC {1,2,3}}. 


Let the set int S* = {x eS: XXX, > 0} and let the set OS” = S* \int S” be the interior 
: 2 : 

and the boundary of the simplex S~, respectively. Let e= (1, 0,0), e,= (0, 1,0), 

e,= (0, 0, 1) be the vertexes of the two-dimensional simplex. 

Theorem 1. For the SCSO W (20), the following assertions true: 


(1) The face Pay ' Dr =e 


; 2 ; ; 
13 ; of the simplex S~ are invariant sets; 
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(ii) Fix(W ) = {e,,e,,e,}; 
(ii) fae [-1, 0) , then @, is an attracting point, €, is a saddle point and @, is a repelling point; 
fae (0,1] , then €, is a repelling point, ©, is a saddle point and €, is an attracting point. 


Proof: (i) Obviously. 


(ii) To find the fixed points we consider the equation W (x) = X, that is the following system of 


equations 


x, = x (I-ax,)(1-$5 J 


x, =x,(1-ax,)(1+ax,), @D 


2 
a 
x, = [14 $5 |(lrax) 


(a) Let X, = (0). Then from the second equation of (21) it follows 
X, =X, —AX,X, => X, =0 and 1 =1-ax,. 
It is clear that if x, = (0 then we have Xx, = 1. If Xx, > 0 then from the last equation one has 


x,=0> x, =1. 


We take solution X, = 1, then it follows that x, = 0. Consequently, if x, = 0), we obtain the fixed 


points €,, e.. 
Similarly, in the case X, = 0 and x, = 0, we have the fixed points €,, €, and €,, €,, respectively. 


(b) Suppose that X,X,xX, F (0). Then from the system (21), one has 
1=(1-ax)(1-$5 } 
2 
1=(1-ax, )(1+ax,), @2) 
a 
1=[14 $5 J(r+ax) 


Ifae [-1, 0) , then from the first equation of (22) we get 
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a a 1 
t= (1=ax(1- $55 | > Bh = tn 


- 
But from QE [-1,0) and from the last equation it follows that xX, = Q and = Q, this 
contradicts to X,X,X, # Q. 


Ifae (0, 1] , then from the third equation of (22) we have 


a a 1 
i= [1s Sx, Jus on 5% = (3 + «) 
But from ad € (0,1] and from the last equation it follows that Xx, = Q and xX, = 0, this contradicts 
to X,X,x, #0. 
Consequently, we have that Fix (W) int S* =O. 


(ii1) To find the type of fixed point of the SCSO (20), we rewrite it in the form 


i =x (I-ar)(1-S 0-5 =) 


. = x, (1+ ax, )(1-a(1-x, —x,)). 


(23) 


where (x,,x,) Ee (x, y) :x,y20,0<x+y< 1 and X,,X, are the first two coordinates of 
a point lying in the simplex S$ a 


a 
The Jacobian of the operator (23) at a fixed point €, has the following eigenvalues A, =l1+-, 


A, =1+4, at the €, has the following eigenvalues A, =l-a, A, =1+4 and at the €, has the 


a 
following eigenvalues A, =|- a, A, =|- 2 Therefore, it follows that: 


Ifae [-1, 0) then, the fixed point €, is a attracting point, €, is a saddle fixed point and the fixed 


point @, is a repelling point. 
If ad =O then, fixed points €,, €, and €, are non-hyperbolic points. 


If ae (0,1] then, the fixed point €, is a repelling point, €, is a saddle fixed point and the fixed 
point €, is a attracting point. 


The theorem is proved. 
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